Phases, transitions, and boundary conditions in a model of interacting bosons 
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We carry out an extensive study of the phase diagrams of the extended Bose Hubbard model, with 
a mean filling of one boson per site, in one dimension by using the density matrix renormalization 
group and show that it contains Superfluid (SF), Mott-insulator (MI), density- wave (DW) and 
Haldane-insulator (HI) phases. We show that the critical exponents and central charges for the 
HI-DW, MI-HI and SF-MI transitions are consistent with those for models in the two-dimensional 
Ising, Gaussian, and Berezinskii-Kosterlitz-Thouless (BKT) universality classes, respectively; and 
we suggest that the SF-HI transition may be more exotic than a simple BKT transition. We show 
explicitly that different boundary conditions lead to different phase diagrams. 
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Over the last two decades, experiments on cold alkali 
atoms in traps have obtained various phases of correlated 
bosons and fermions [l[. Microscopic interaction param- 
eters can be tuned in such experiments, which provide, 
therefore, excellent laboratories for studies of quantum 
phase transitions, such as those between superfluid (SF) 
and Mott-insulator (MI) phases in a system of interact- 
ing bosons in an optical lattice 0-0] ■ The SF and MI 
phases are, respectively, the protypical examples of gap- 
less and gapped phases in such systems; in addition, it 
may be possible to obtain exotic phases, e-g^ density- 
wave (DW) 041, Haldane-insulator (HI) 0, lid ], and su- 
persolid (SS) [11| phases, in systems with a dipolar con- 
densate of 52 Cr atoms [l2|. The first step in developing 
an understanding of such experiments is to study lat- 
tice models of bosons with long-range interactions fl3j |: 
the simplest model that goes beyond the Bose-Hubbard 
model with onsite, repulsive interactions 0, Q is the ex- 
tended Bose-Hubbard model (EBHM), which allows for 
repulsive interactions between bosons on a site and also 
on nearest-neighbor sites [3, H]. The one-dimensional 
(ID) EBHM is particularly interesting because (a) quan- 
tum fluctuations are strong enough to replace long-range 
SF order by a quasi-long-range SF, with a power-law de- 
cay of order-parameter correlations @, and (b) boundary 
conditions play a remarkable role here, in so far as they 
can modify the bulk phase diagram [3, [13, EH EH > m a 
way that is not normally possible in the thermodynamic 
limit 0. 

We present extensive density-matrix-renormalization- 
group (DMRG) [3| studies of the phases, transitions, and 
the role of boundary conditions in the ID EBHM with 
a mean filling of one boson per site. Although a few 



DMRG studies @, i, E3, EI El have been carried out 
earlier, none of them has elucidated completely clearly 
the universality classes of the continuous transitions in 
the ID EBHM; nor have they compared in detail, the 
phase diagrams of this model with different boundary 



conditions. Our study, which has been designed to obtain 
the universality classes of these transitions and to exam- 
ine the role of boundary conditions in stabilizing differ- 
ent phases, yields a variety of interesting results that we 
summarize qualitatively below: If the filling p = 1 , in the 
thermodynamic limit, we obtain three types of phase dia- 
grams, PI (Fig. Ufa)), P2 (Fig.rjjb)), and P3 (Fig. [17c)), 
for boundary conditions BC1 (number of bosons N = L, 
where L is the number of sites), BC2 (N = L + 1), and 
BC3 (N = L and additional chemical potentials at the 
boundary sites, as we describe below), respectively. For 
PI we obtain SF, MI, and DW phases and we confirm, 
as noted earlier [6j, that the SF-MI transition is in the 
Berezinskii-Kosterlitz-Thoulcss (BKT) [I?} universality 
class and the SF-DW has both BKT and two-dimensional 
(2D) Ising characteristics; at sufficiently large values of 
the repulsive parameters, the MI-DW transition is of first 
order @. The phase diagram P2 displays SF, HI, and 
DW phases, with continuous SF-HI and HI-DW transi- 
tions in BKT and 2D Ising universality classes, respec- 
tively. The phase diagram P3 has SF, MI, HI, and DW 
phases, with SF-MI, MI-HI, HI-DW in BKT, 2D Gaus- 
sian (l8j, and 2D Ising universality classes, respectively; 
the SF-HI transition is more exotic than a simple BKT 
transition; at sufficiently large values of the repulsive pa- 
rameters, there is a first-order MI-DW phase boundary. 
The ID EBHM is defined by the Hamiltonian 

n = -t^[ala^ 1 +h.c)+—^2hi(h i -l)+V^2 



n-im+i, 



(1) 

where t is the amplitude for a boson to hop from site i 
to its nearest-neighbor sites, h.c. denotes the Hermitian 
conjugate, at, Oj, and hi = a\ai are, respectively, boson 
creation, annihilation, and number operators at the site 
i, the repulsive potential between bosons on the same site 
is U, and V is the repulsive interaction between bosons 
on nearest-neighbor sites. We set the scale of energies by 
choosing t = 1. We work with a fixed value of the filling 
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p and focus on p = 1, i.e., a filling of one boson per site 
in the thermodynamic limit, which we realize in three 
different ways, to highlight the role of boundary condi- 
tions, by using the three different boundary conditions 
BC1, BC2, and BC3. In all these three cases we have 
open boundaries; in BC3, we include boundary chemical 
potentials pi = —It and p r = 2t at the left and right 
boundaries, which are denoted by the subscripts I and r, 
respectively. We study this model by using the DMRG 
method that has been described in Ref. [(| . To character- 
ize the phases and transitions in this model, we calculate 
the following thermodynamic and correlation functions 
and the entanglement entropy (familiar from quantum 
information): 

R dw {\i-j |) = (-l^iSfkShi), 
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where Shi = hi — p, the correlation functions for the 
DW, HI and SF phases are R dw {\i - R s tring(\i - 
RsF(\i — j\)> respectively, is the density-density struc- 
ture factor at wave number k = tt, n(k = 0) the mo- 
mentum distribution of the bosons at k — 0, G^ G and 
G£ G are neutral and charge gaps 0, El(N) and E\{N) 
are the ground-state and first-excited-state energies, re- 
spectively, for our system with L sites and N bosons, 
£l is a system-size-dependent correlation length; Sx(Z) 
is the von-Neumann block entanglement entropy, A,; are 
the eigenvalues of the reduced density matrix for the 
right block, of length I, in our DMRG @, and n sta tes 
is the number of states that we retain for this de nsity 
matrix; note that Sl{1) = cA at a critical point [20(, 
with c the central charge and the log-conformal distance 
A = (log [2^ sin(^)])/6; a and B depend on the criti- 
cal value V c , at which the concerned transition occurs; 
V c depends on U. The order parameters for the DW 
and HI phases are Od w = yTimuZjjZ^ 



Rdw{\ i — and 



for the DW, HI and SF phases, we also use the Fourier 
transforms of the correlation functions, which we denote 
by Rdw(k), Rstring(k), J?si?(fc), respectively. Our cal- 
culations have been performed with 100 < L < 300, a 
maximal number n max = 6 of bosons at a site, and up 
to n sta tes = 256 states in our density matrices; we have 
checked, in representative cases, that, for the ranges of U 
and V we cover, our results are not affected significantly 
if we use n max = 4 and n states = 128. 

In Figs. Q] (a), (b), and (c) we depict our DMRG phase 
diagrams, in the (U, V) plane, for the ID EBHM with 
BC1, BC2, and BC3 boundary conditions, respectively. 
These phase diagrams show MI (ochre), SF (purple), HI 
(red), and DW (green) phases and the phase boundaries 
between them; the points inside these phases indicate 
the values of U and V for which we have carried out 
our DMRG calculations. In the ranges of U and V that 
we consider here, all transitions are continuous; at larger 
values of U and V than those shown in Fig. [TJ the MI- 
DW transitions become first-order, as shown for BC1 in 
Ref. Q; the direct MI-DW transition, at large values of 
U and V, is also of first order with the BC3 boundary 
condition. Note that MI (HI) phases appear with bound- 
ary conditions BC1 and BC3 (BC2 and BC3). The MI 
phase occurs at a strict commensurate filling (BC1 with 
N = L) and is characterized by a charge gap, which 
arises because of the onsite interaction U that penalizes 
the presence of more than one boson at a site. With the 
boundary condition BC2, iV = L+l, so there is one more 
boson in the system than in case BC1; this extra boson 
can hop over the commensurate, BC1, MI state, produce 
gapless excitations in the upper Mott-Hubbard band, and 
thereby destroy the MI phase, which is replaced by SF 
or HI phases. 

The origin of the HI phase in the ID EBHM can be 
understood qualitatively as follows: If p = 1 and U ^> t, 
the MI phase has one boson at most sites because the 
large value of U suppresses boson-number fluctuations; 
as U begins to decrease, these fluctuations are enhanced 
and, to leading order in |i|/J7, arise from holes at some 
sites and pairs of bosons at other sites (occupancies of 
more than 2 bosons at a site are energetically disfavored 
unless U decreases). If we restrict ourselves to 0, 1, and 2 
bosons at every site, the ID EBHM can be mapped onto a 
spin-1 chain, in which the three values of spin projection 
correspond to the occupancies 0, 1, and 2. This spin-1 
model, an antiferromagnetic chain, with full SU (2) sym- 
metry broken by the extended interaction V in Eq.([T]), 
has a Haldanephase, whose analog in the ID EBHM is 
the HI phase [9(; this phase also has gapless edge modes, 
which destroy the charge gap at strict commensurate fill- 
ing (BC1). However, with an extra boson (BC2), these 
gapless modes get quenched and so we obtain an HI phase 
with a gap. 

By applying the boundary chemical potentials pi and 
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FIG. 1. (Color online) Phase diagrams in the {U,V) plane 
for the ID EBHM with the following boundary conditions (a) 
BC1 (N = L), (b) BC2 (N = L+l), and (c) BC3 (N = L) and 
fir = — fii = 2 showing Mott-insulator (MI ochre), superfluid 
(SF purple), Haldane-insulator (HI red), and density- wave 
(DW green) phases and the phase boundaries between them; 
in this range of U and V all transitions are continuous; at 
larger values of U and V the MI-DW and HI-DW transitions 
become first-order. 



fi r in BC3, we introduce a gap in the edge states in the 
HI; however, because the system is exactly at a com- 
mensurate filling N = L, the MI phase is also present. 
Thus, we obtain both the MI and the HI phases with the 
BC3 boundary condition. Note also that the HI phase 
exists because of boson-number fluctuations about the 
Mott state; so we expect the HI to give way to the MI 
phase at large values of U at which such number fluc- 
tuations are suppressed, in much the same way as the 
SF yields to the MI at large values of U. In the large- U 
and large-]/ regions, which we do not show in the phase 
diagrams of Fig. [TJ there is a direct, first-order, MI-DW 
transition Q. 

We now obtain the universality classes of the contin- 
uous transitions in the phase diagrams of Fig. Q] It has 
been noted in Ref. @ that the MI-SF transition in Fig. Q] 
(a), with BC1 boundary conditions, is in the Bcrezinskii- 
Kosterlitz-Thouless (BKT) class, whereas the SF-DW 
has both BKT and two-dimensional (2D) characters, in 
as much as the SF correlation length shows a BKT diver- 
gence but the DW order parameter decays to zero with 
a 2D-Ising, order-parameter exponent; this transition 
should be in the universality class of the Wess-Zumino- 
Witten SU(2)i theory like the critical point that lies be- 
tween the BKT phase and the antiferromagnetic phase 
in the spin— 1/2, XXZ antiferromagnetic chain [2l|. The 
MI-SF and SF-DW phase boundaries in Fig.[T](a) merge, 
most probably at a bicritical point, beyond which the MI- 
DW transition is first order |6j. We concentrate here on 
the HI-SF, HI-DW, HI-MI transitions, all of which can 
be found in Fig.[T](c), with the boundary condition BC3, 
for which we present results in Fig. [2] we have checked 
explicitly that the universality classes of the HI-SF and 
HI-DW transitions are the same with both BC2 and BC3 
boundary conditions. 

In Fig. [2] (eti) we plot the scaled DW order param- 
eter OdwL 13 /" versus V near the representative point, 
V c ~ 3.86 and U = 6, on the HI-DW phase boundary 



for L = 200 (red curve), L = 150 (blue curve), and 
L = 100 (green curve); the inset gives a plot of Od w L^/ v 
versus (V — V c ) v L\ Fig. [2] (02) presents a similar plot and 
inset for the scaled structure factor S{ji)L 2 ^l v at wave 
number k = tt; in Fig. [5] (013) we plot the von-Neumann 
block entanglement entropy Sl(1) versus the logarithmic 
conformal distance A for L = 220 and a range of values 
of V in the vicinity of VJJJ = 6) = 3.86. The values 
of the exponents j3 and v and the value of the central 
charge c that we obtain from Figs. [2] (ai) — (03), namely, 
/3 = 1/8, v = 1, and c = 0.52, are consistent with their 
values in 2D Ising model. If we restrict the occupancy to 
one boson per site, in the limit U, V — s- 00 with fixed U/V, 
the DW state can be represented as | . . . 20202020 . . . ) in 
the site-occupancy basis; this state is doubly degenerate 
and clearly breaks the translational symmetry of the lat- 
tice by doubling the unit cell; in contrast, the HI state 
has the translational invariance of the lattice. Given this 
double degeneracy of the DW state, we expect that, if 
the HI-DW transition is continuous, then it should be in 
the 2D Ising universality class; as we have shown above, 
our DMRG results for the HI-DW transition in the ID 
EBHM arc consistent with this expectation. Note that 
the string order parameter O s t r ing is nonzero in both 
the HI and DW phases; thus, it is not the order pa- 
rameter that is required for identifying the universality 
class of the HI-DW transition. We find, furthermore, the 
string-order-parameter correlation length is so large in 
the HI phase that, given the values of L in our study, it 
is not possible to get a reliable estimate for this correla- 
tion length. 

To investigate the critical behavior of the MI-HI tran- 
sition, we plot R s tring(k = 0)L 2 ^ versus V, in Fig. [2] 
(61), whose inset shows a plot of R s tring(k = 0)I/ 2/3 / l/ 
versus (V — V c ) v ' L. These scaling plots indicate that, for 
[7 = 6, the MI-HI critical point is at V c ~ 3.5 with ex- 
ponents j3 — 0.35 and v = 1.18; this value of j3 is consis- 
tent with the Gaussian-model [18[ (superscript G) result 
fi G = l/v8. Our value for v follows from the plot of 
L/S.L versus V in Fig. [2] (62), whose inset shows a plot of 
L/£i versus (V — V^L 1 ^ . For the central charge we ob- 
tain c = 1.02, which is consistent with c = 1 given our 
error bars, from the plot of Sl(1) versus A for L = 220 in 
Fig. [U (63). As we have mentioned above, for large U we 
can restrict the occupancies in the ID EBHM to 0, 1, and 
2 bosons per site to obtain a spin-1 model. The analog of 
the MI-HI transition in this spin-1 model has been shown 
to be in the Gaussian universality class [l8| , which has a 
fixed exponent /3 = l/\/8 and central charge c = 1, but 
an exponent v that changes continuously along the phase 
boundary. We find indeed, from calculations like those 
presented in Figs. [2] (61) — (63), that j3 and c do not change 
as we move along our MI-HI phase boundary in Fig. [1] 
(c), but v does; we obtain, in particular, v = 1.18, 1.22, 
and 1.36 for U = 6, 5, and 4, respectively. 

The critical behavior of the SF-HI transition follows 
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from the plot of R str m g {k = Q)L 2 P/ V versus V, in Fig. [2] 
(ci), and the the plot of L/£i versus V, in Fig. [2] (c 2 ), 
whose inset shows a plot of F versus D (see Eq.[5]). These 
scaling plots indicate that, for V = 0.9U, the SF-HI crit- 
ical point is at V c ~ 1.7, where the correlation length di- 
verges with an essential singularity as it does at a BKT 
transition [l7j |. Thus, we cannot define the exponent v; 
however, we can define the exponent ratio fi/v, which 
governs how rapidly the string order parameter vanishes 
as we approach the SF-HI phase boundary from the HI 
phase; for this ratio we obtain the estimate f3/v ~ 0.9. 
We find the central charge c = 0.94, which is consistent 
with c = 1 given our error bars, from the plot of Si, (I) 
versus A for L = 280 in Fig. [2] (c 3 ). The field theory for 
this SF-HI transition is likely to be non-standard because 
the SF phase has algebraic U(l) order, whereas the HI 
phase has a non-local, string order parameter. An exam- 
ple of a similar transition is found in the phase diagram 
of the bilinear-biquadratic spin-1 chain, in which the Hal- 
dane phase (the analog of the HI phase here) undergoes 
a transition to a critical phase, as we change the ratio 
of the bilinear and biquadratic couplings; this transition 
is known to be described by an SU(3)± Wess-Zumino- 



Wittcn theory [22j. The gaplcss phase in our model is 
quite different from that in the spin-1 chain; and there 
is no SU(3) symmetry at the SF-HI transition in the ID 
EBHM. Nevertheless, it is likely that the SF-HI phase 
transition in our model is of a non-standard type, which 
is similar to, but not the same as, the SU(3)\ theory 
mentioned above; this point requires more detailed inves- 
tigations that lie beyond the scope of our paper. Note, 
furthermore, that the SF-HI transition here is not fine- 
tuned, unlike the one above in the spin-1 chain with bi- 
linear and biquadratic couplings, in so far as the SF-HI 
transition occurs along an entire boundary in the phase 
diagram (Fig. Q] (c)). 

For the sake of completeness we also show, for the BC3 
boundary condition, the BKT nature of the SF-MI transi- 
tion in Fig.Q](c), as noted for the BC1 case in Ref. 0, by 
presenting plots of the scaled charge gap F versus D (see 
Eq. [2]), the k = value of n(k = 0), and the block entan- 
glement entropy Sl(1) versus the log-conformal distance 
A in Figs.[2](<ii), (d 2 ), and (d 3 ), respectively. These plots 
demonstrate that, at the MI-SF transition, the correla- 
tion length (or inverse charge gap) has a BKT-type essen- 
tial singularity (Figs. [2] di), the SF-correlation- function 
exponent r) = 1/4, and the central charge is c = 0.96, 
which is consistent with c BKT = 1 given our error bars. 
It is difficult to give precise error bars for the exponents 
and central charges that we have calculated because there 
are systematic errors, which include those associated with 
the finite values of n rnax and n sta tes these errors are hard 
to estimate. 

At lower values of U and V than those shown in Fig. [IJ 
it has been suggested that the ID EBHM can show a 
supersolid (SS) phase between SF and DW phases (see, 



e.g., Ref. 14 1). However, it is not easy to distinguish such 
an SS phase from a region of first-order phase coexistence 
between SF and DW phases in a DMRG calculation, with 
a fixed number of bosons; we discuss this elsewhere [23| . 

We have performed the most extensive DMRG study 
of the ID EBHM attempted so far. Our work, which ex- 
tends earlier studies [1, 0, 0, E3] significantly, has been 
designed to investigate the effects of boundary conditions 
on the phase diagram of this model, to elucidate the na- 
tures of the phase transitions here, and to identify the 
universality classes of the continuous transitions. Our 
study shows that the HI-DW, MI-HI, and SF-MI tran- 
sitions are, respectively, in the 2D Ising, Gaussian, and 
BKT universality classes; however, the SF-HI transition 
seems to be more exotic than a BKT transition in so far 
as the string order parameter also vanishes where the SF 
phase begins. We find that the different boundary condi- 
tions BC1, BC2, and BC3, which lead to the same value 
of p in the thermodynamic limit, yield different phase 
diagrams; this boundary-condition dependence deserves 
more attention than it has received so far, especially from 
the point of view of the existence of thermodynamic lim- 
its [H. We hope our work will stimulate more theoretical 
and experimental studies of the ID EBHM and experi- 
mental realizations thereof. 
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